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We present a consistent second order perturbation theory for the lowest-lying condensed modes of 
very small, weakly-interacting Bose-Einstein condensates in terms of bare particle eigenstates in a 
harmonic trap. After presenting our general approach, we focus on explicit expressions for a simple 
three-level system, mainly in order to discuss the analogy of a single condensate occupying two modes 
of a trap with the semi-classical theory for two-mode photon lasers. A subsequent renormalization 
of the single-particle energies to include the dressing imposed by mean fields demonstrates clearly 
the consistency of our treatment with other kinetic approaches, 
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I. INTRODUCTION 

The kinetic theories of Bose-Einstein condensation in a trap can be divided into two general categories. The fully 
quantum approaches have been based either on a perturbative master equation treatment for the many-body density 
matrix [0 , or on a single Fokker-Planck equation for the nonequilibrium dynamics of the entire system . Alternative 
approaches to the nonequilibrium dynamics of trapped Bose-condensed gases are essentially of a perturbative nature, 
based on a suitably truncated coupled equation of motion hierarchy for normal and anomalous averages . The 
growth of condensation was first studied numerically by Gardiner and co-workers [p^ , pT| , based on the theory developed 
in 01, and their findings were in qualitative agreement with the description of Kagan, Svistunov and Shlyapnikov 
|12|. Independent studies by Stoof and co-workers, based on similar approximations, have produced growth curves 
1 13[ which are in very good agreement with those of Gardiner et al. . 

The large condensates typically produced in experiments lead to large mean field potentials which significantly dress 
the single-particle eigenenergies of the trap potential, thus making it convenient to work in a suitably dressed basis. 
The conventional theoretical picture adopted is thus that of a very large single-mode condensate interacting with a 
large number of quasiparticles and higher-lying excited atoms, whereas it is not a priori necessary that single-mode 
condensation will correspond to all experimental conditions Instead of working in the usual 'condensed matter' 

approach of a condensate and a set of quasiparticle excitations which are, by definition, orthogonal to it (which is 
suitable for large condensates) we focus our description in terms of a more 'quantum-optical' approach, i.e. in terms 
of individual modes of the trapping potential |^|^ , for which one must in principle consider the condensate spanning 
a large number of modes, with the same modes being simultaneously occupied by non-condensate atoms. Within 
such a picture, this paper explicitly discusses the simplest deviation from single-mode condensation in the context 
of a bare single-particle basis, in a manner analogous to existing treatments. Formulating the problem in terms of 
bare particle eigenenergies can only be useful in the case of extremely small, dilute, weakly-interacting condensates, 
when the mean field effects are small enough that their induced shifts of the single-particle energies can be treated as 
perturbations. Such an approach is nevertheless beneficial for discussing the inherent multi-mode nature of trapped 
condensates and their relation to the corresponding (semi-classical) theory of multi-mode photon lasers. 

The first part of this paper (Sec. II) reviews in the usual manner the formal development of the coupled equation 
of motion approach. To simplify the physical picture and bring out the underlying structure in a clear manner, we 
focus on a very simple system consisting of the three lowest trap eigenstates, for which we derive explicit equations 
of motion in the Popov approximation in Sec. Ill (with corresponding off-equilibrium contributions given in 
Appendix A, and some further clarifications in Appendix B). By explicitly discussing the interplay of two coupled 
condensed modes, we show how our treatment reduces to the Hartree-Fock theory for binary condensates [^6| (Sec. 
IV Al). An important advantage of formulating the theory in terms of bare single-particle eigenstates is that it allows 
us to discuss in detail the analogy between our equations for the coupled dynamics of two condensed modes, and the 
corresponding ones arising in the semi-classical treatments of two-mode photon lasers (Sec. IV A2). We show that, 
the extent to which this analogy can be drawn for the inherent multi-mode nature of a single condensate, depends 
critically on the choice of the single-particle basis (i.e. bare or dressed by various types of mean fields) in terms of 
which the analysis is carried out. Our discussion here is distinct from conventional analogies based on two different 
condensates which are spatially separated in different spin states p8[, or outcoupled by radiation applied at 
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two different frequencies [g^. Sec. IV B compares our approach to conventional kinetic treatments, where we show 
exphcitly that our theory reduces to the multi-mode kinetic treatment of Walser et al., upon shifting our description 
to a basis in which the single-particle eigenenergies become dressed by the usual Hartree-Fock-Bogoliubov (HFB) 
mean field potentials. This suggests that, contrary to an implication in such dynamic equations (and coUisional 
integrals) are inevitably basis-dependent. 



II. THE COUPLED EQUATION OF MOTION APPROACH 

Consider a sufficiently dilute, weakly-interacting partially Bose-condensed trapped gas with a binary-interaction 
hamiltonian 

H = ^ Sf/'''=aJ;as + V'rsmnalalamO.n (1) 

rs rsmn 

Here S''"''^ — — (?i^V^)/(2m) + Vtrap(r) contains both kinetic energy and trapping potential and Vrsmn rep- 
resents the symmetrized form of the interaction potential between a pair of particles, defined by Vrsmn = 
i |(rs|y|mn) -|- (rs|V^|nm)|, where \i) = ipi{r) denotes a single-particle eigenstate of the trap. The single-particle 

operators are related to the Bose field operator >I'(r, t) via \l/(r, t) = J^i V'i(r)ai(0. We assume the system to be in 
a symmetry-broken phase and hence express the single-particle operators ai{t) as pG] 

a^{t) = {a,) + {a, - (a,)) = z^{t) + c.,{t) (2) 

This allows, in general, for a coherent mean-field amplitude Zi to form in a number of low-lying trap levels. We can now 
formulate a non-equilibrium theory for the coupled evolution of condensate mean field amplitudes and fluctuations 
about these values, based on a suitably truncated hierarchy of coupled equations of motion and appropriate decoupling 
approximations. 

Possibly the most direct approach for studying the dynamics in a closed system is based on solving an appropriate 
set of such equations self-consistently, in terms of exact interatomic potentials; such an approach has been discussed, 
for example, in In a realistic system, the number of trap eigenstates will be very large, making such a procedure 
computationally very demanding. Since the most interesting dynamics take place in the low-lying levels, in this paper 
we have chosen to restrict our analysis to such levels (although such treatment will also implicitly yield the behavior 
of high- lying thermal levels). Accepting a distinction between low- and high- lying levels allows us to adiabatically 
eliminate all high-lying levels appearing as intermediate states in the equations for the evolution of averages of low- 
lying states; this procedure is known to lead to the renormalization of the exact (single-vertex) interatomic potential 
to an effective two-body one (over high- lying states), as discuscd in [^. We thus arrive at the situation where the 
hamiltonian of the system still has the general form of Eq. (1), but with the single- vertex interatomic potential V 
replaced by an effective two-body T-matrix, T, over high-lying states, with the simultaneous restriction of all bare 
trap eigenstates being summed over low-lying levels. For sufficiently dilute systems at low temperature, this restricted 
effective two-body interaction is approximately equal to the full two-body T-matrix, T^^, giving the scattering of 
two particles in vacuum. We thus approximate T in terms of T^^, by ensuring that purely two-body effects due 
to collisions occuring in vacuum are not double-counted (for a detailed discussion of the relation of these effective 
interactions and the renormalization required to avoid double-counting, the reader is referred to | |2l| , ^ ). We note 
that it is precisely the quantity T^^ which corresponds, in three dimensions, to the usual binary s-wave scattering 
length pseudopotential. 

We are interested in working out in a self-consistent manner the evolution of condensate population and incoherent 
ffuctuations about this value, and for consistency we work throughout this paper with averages of two single-particle 
operators. The Heisenberg equation of motion for such a general product of two operators is given by (setting h = 1) 



d 



r 

+ [Tjsmrialdlamdr) - T^*^™,. (aL^a^aj) } (3) 

rms 

where the primes indicate summation over low-lying levels. Consistent application of second order perturbation theory 
in the weakly-interacting limit (i.e. when the system can be well described in terms of single-particle wavef unctions) 
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should yield correct expressions for energy level shifts and population damping. To proceed with our treatment, we 
must determine whether there are any quantities (fluctuations) which evolve faster than others (mean fields), so that 
the former can be adiabatically eliminated from the first order expressions. For example, in the usual rate equation 
treatments, one eliminates all off-diagonal normal and anomalous averages in favour of (diagonal) populations; one 
often speaks of coherences damping out faster than populations due to the coupling of the system to its environment. 
This gives rise to a set of equations to second order in the effective potentials, coupling populations to populations, 
in what is often termed the secular approximation p^ . 

One could, however, argue that the choice of which low-lying averages can be adiabatically eliminated depends 
on the basis employed for the description of the system, i.e. essentially on whether the mean field energy shifts are 
correctly taken into account or not. Starting from bare trap eigenenergies, one would not expect any normal or 
anomalous averages to be slowly-evolving; hence, the correct second order expression of Eq. (3) in a bare particle 
basis can be obtained by adiabatic elimination of the entire quantity {a^a'^aa) by means of its respective equation of 
motion. This is the procedure adopted in the main part of this work, which gives rise to complex equations appearing 
to contain additional terms when compared directly to similar treatments. However, a subsequent transformation to 
a dressed single-particle basis (i.e. dressed eigenenergies) shows clearly that such terms drop out from the respective 
equations of motion for populations in dressed eigenstates, as anticipated. Nonetheless, an approach in terms of bare 
single-particle eigenenergies allows one to draw important analogies between multi-mode condensation and multi-mode 
laser theory. Proceeding thus with the treatment in a bare basis, and taking the operator S''"'"'^ to be diagonal, we 
obtain 

i-^ialalamUr) = {uj,n + ujs){a\alar,iar) + Fi{T, a, a'';t) (4) 

where tUi correspond to bare trap energies and Fi{T, a, a^; t) defines the collisional evolution of (a|aja„iar) in such a 
basis by Q 

Fi{T, a, a^) = ^ TrarH{a\alaiat) - ^ Tpqis{a\,a\arnar) 

It pq 

+ Tprit {a\,a\alamaiat) + Tprait{a\,a\alaraiat) 
pit pit 

- ^ Tpqis {al,a\a\aiarnar) - ^ Tpqu {ala\alaiarnar) (5) 

pql pql 

Assuming real eigenvalues, we obtain the following exact integral relation 
d 



rm.t •'*0 



(6) 



where e.c. stands for the exchange conjugate (i.e. conjugate expression with labels i and j interchanged). By using 
the definition a, — Zi + Ci, we obtain all second-order collisional terms of our approach. However, for these to be 
useful, we must express them in terms of a closed system of equations, by imposing suitable approximations. Firstly 
we decouple averages containing more than two single-particle operators via 

(clclcmCn) ~ {clcm) {clc„) + {clcn) {c\cm) + {clc\) {CrnCn) (7) 



and 



{clclclcqClCt) ~ {clcq) ((cJq) (c^Ct) + {clct) {clci)) 

+ {clci) {{clcq){clct) + {clct){clcq)) + (c],Ct) {{clcq) {c\ci) + {clci) {0%)) 

+ (cqQ) ((ctct)(ctct) + {clcl){clct) + (ctct)(ctct)) 

+ {CqCt) ((ctct)(ctQ) + (ctct)(ctQ) + (cj ) (ct C, ) ) 

+ {cict) {{clcl){clcq) + {clc\){clcq) + {clc\){clcq)) (8) 
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Based on our formulation in terms of effective interactions, we further impose the Markov approximation which takes 
the quantity Fi{T,a,a^) out of the integrand and defines the intermediate propagators. This imphes the assumption 

that the operators Oj evolve freely between collisions via di{t ) — e+*"»(*~* 'iai{t). Thus, the second order contributions 
of Eq. m) acquire the general form 



T2 



E 



T,, 



^^g-»(A^...)(t-t') j T...Fl..{t) 



+ e.c. 



(9) 



where the above quantity T...Fi...[t) corresponds to Fi defined by Eq. (5), and this notation has been used to indicate 
that the dotted indices of the second T are the same indices as the ones appearing in the exponential of the integrand, 
as a result of the Markov approximation. Following the notation of Gardiner et al. 0], we now rc-write the above 
integral as 



/>OC 

lim / dre^^^'^^^"')^ = 5^p\I^w) = 7r^(Aa;) ± iV 

17^0+ Jo 
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(10) 



where V (■^) corresponds to the principal value of the integral and the upper limit of integration has been approxi- 
mated by (t — to) oo, based on the usual assumption that successive coUisional events are well separated in time. 
Having discussed our approximations, we can now write down the coupled rate equations for coherent, incoherent and 
total populations. 



III. RATE EQUATIONS FOR MULTI-LEVEL CONDENSATION 



In this section, we focus on the application of the above methodology to a simple three-level system which we 
discuss within the Popov approximation |]l5| , in which anomalous averages of the non-condensate are ignored in 
the final expressions. We are well aware that this system is rather idealized and do not claim that it will accurately 
reproduce the entire dynamics of condensed and thermal atoms. In fact, corrections beyond Popov may be significant, 
as shown in this context by Walser et al. Q (see also [p5[-p7t), and we defer their explicit discussion to a subsequent 
paper (but see also Sec. IV. B). In the current paper we are mainly concerned with addressing the simplest departure 
from single-level condensation, and how this might affect the dynamics of low-lying levels. To this aim, such a system 
is ideal for a simple comparison with other kinetic theories ]l|J^,^,[7|,p|jl^] . Perhaps more significantly, such a small 
system will enable us to discuss in an explicit manner the desired analogy of the multi-mode nature of a single 
condensate with multi-mode semiclassical laser theory Our notation is as follows: The total population Ni in 

level i is written as a sum of coherent and incoherent rii = pa populations, via 

TV, = {ala^) = (c|c,) + \z,f = n, + \z^f (11) 

with off-diagonal coherent and incoherent matrix elements respectively defined by — {z*Zij and pij — {c^jCi). 



A. Evolution of Total Populations 

We start with expressions for the coupled evolution of total populations in the three lowest-lying bare levels under 
the assumption that all of them may exhibit partial condensation. Since we are dealing with a closed system, the 
evolution of populations in these levels trivially satisfy 



dNa dN2 1 fdNi 



dt dt 2 V dt 

where 



(12) 



dt 

Here we have defined the collisional redistribution rate 



I [1^0211 (PIO + Clo) P12 - C.c] + { (i?T|„,,,|. + Q'l^ll) + C.C.} (13) 
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where 



= 2ro2ii [{No + 1){N2 + 1)N^ - NoN2{Ni + if 
-ro2ii(l + A^o + N2)\zi\^ + Rt.,,,,., 



ro211 = IT0211P lim / d^e±^(-o+-2-2^i±»r,)r ^ jTo^n (c^q + ^2 " 2c^i 



(14) 



(15) 



and the 'tilde' in ^T^^^^p and Qifjl denotes contributions which are off-diagonal in p or whose explicit expressions 
can be found in Appendix A. In the context of the well-known Boltzmann scattering factors, the appearance of the 
additional term ~ ro2ii|zi|^ may appear somewhat perplexing. We stress that such contributions only arise because 
of our choice to formulate our description in terms of a bare basis, and we explicitly show in Sec. IV B that such terms 
'disappear' in the usual formulation in terms of a self-consistently dressed basis (where they become incorporated in 
the dressing of the bare trap eigenstates). 

Unlike the evolution of total populations, the dynamics of condensed and thermal atoms will also depend on 
processes of the general form TijjiTkuk which will be intrinsically dispersive. For simplicity, we shall henceforth limit 
our discussion to the case of only two partially condensed levels, namely 2:0,21 7^ 0, with level 2 being the lowest 
purely incoherent level treated here. 



B. Evolution of Condensate Populations 



The method used to obtain the evolution of condensate populations, is essentially the same as the one used for total 
populations (Eq. (|^)).For ease of our subsequent comparison with the semi-classical laser theory, in this section we 
directly obtain equations for the modulus jzip of the condensate mean field amplitude Zi and not for Zi itself. These 
are given by 



dt 



(16) 



where the quantity analogous to Fi of Eq. (5) is obtained from the equation of motion of the entire quantity z* (aJdmOr) 
and not just {d\dmCLr) (see Appendix B). Thus, to second order in the effective interaction, we obtain 



d\zi 



dt 



C.C. 



C.C. 



(17) 



with all coUisional contributions summarized in XulT^), which can be written as 



(18) 



i?™" corresponds to redistribution collisions affecting the level i condensate. Si corresponds to self-interactions within 
level i, and Ei gives the exchange collisions of level i with its nearby levels, whereas the term XuiT'^) corresponds to 
all off-diagonal contributions. For two condensed levels, in our simplified three-level system, the respective expressions 
are 



2|2o|' [{n2 + l)ni-n2{ni + lf] 
i?™" = ro2ii { +4|zorki|' [("2 + l)ni ~ n2(ni + 1)] 
+ koPki|MK + l)-n2] 



(19) 



So = -X 



0000 



ZkoP [(no-f l)2no-n2(no + l)] 
+ \zo\' [(no + I f - nl] 



(20) 



Eq — — 4Aoi 



koP [("0 + l)(ni + l)n-i - noni(ni -|- 1)] 



+ U'\zi\' 



[{no + l)(ni + 1) - noni] 
-I- [(no + l)ni - no(ni + 1)] 



-4Ao22okoP [(no + l)(n2 + l)n2 - non2(n2 + 1)] 



(21) 
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and 



2|zi|2 [(m + 1)710^2 - ni{no + l)(n2 + 1)] ^ 
pcon_r,p ; +2|zo|2|zi|2 [(m + l)n2 - ni(n2 + 1)] 

-\zo\^\zi\H{n2 + l)-n2] 



(22) 



Si = -Xi 



(23) 



£;i = -4A, 



Olio 



kip [(no + l)(ni + l)no - noni{no + 1)] 

, I, 121, 12/ [('^o + l)(ni + l)-no"i] 1 
+ l^il \ +[(ni + l)no-ni(no + l)] / 



-4Ai22i|2;iP [(n-i + l)(n2 + 1)^2 - nin2(n2 + 1)] 



where 



Ayji = \T,j^i\ lim / dre" 



while the expressions for the off-diagonal contributions XaiT'^) can be found in Appendix A. 



(24) 



(25) 



C. Evolution of Non-condensate Populations 

The evolution of the normal averages of levels and 1 is governed by the equation 
d 



dt 



c]ci) = - 



iTi, 



+ c.c. 



(26) 



where Y,,{T^) = Rf - - E, + Y,^{T^) and Rf are defined by 

2 [(no + l){n2 + l)n? - non2(ni + 1)^] 
Rf = ro2ii { +4|zip [(no + l)(n2 + f )ni - non2(ni + 1)] 
+ |2i|M('^o + 1)K + 1) - non2] 



(27) 



[(ni + lfnon2 - n?(no + l)(n2 + f)] 
nth_.p ) +|zo|2 [(ni + l)2n2 - n2(n2 + f)] 
^ " +|zi|2[(ni + l)non2-ni(no + f)(n2 + l)] 

+ |zohzi|2 [(m + i)n2 - ni(n2 + f)] 

The evolution of the non-condensate component of level 2 is given by Eqs. (12)-(14), in the limit Z2 
respective expressions for the off-diagonal contributions Yii{T'^) can be found in Appendix A. 



(28) 



0. The 



IV. INTERPRETATION OF RATE EQUATIONS 



We now turn our attention to the interpretation of the above equations, focusing in Sec. A only on the condensate 
evolution. Ignoring, at first, re-distributional dynamics and off-diagonal contributions, we show how they reduce to the 
Hartree-Fock coupled nonlinear Schrodinger equations conventionally used for describing two-component condensates 
(Sec. IV Al). By additionally considering the redistribution dynamics, we further discuss the analogy of a single 
inherently multi-mode condensate to the semi-classical photon laser theory [ p8[ , with our discussion being limited to 
only two modes for clarity (Section IV A2). By further considering the evolution of uncondensed, and thus also of 
total populations of each level, we compare and contrast our approach to the kinetic theory of Walser et al. which 
is the formal multi-mode kinetic theory closest to the formalism of this paper (Sec. IV B). 
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A. Coupled Condensate Evolution 



1. Hartree-Fock Theory for Binary Condensates 



Let US initially forget about the presence of level 2 in our system and focus on the interactions between the two 
lowest levels in the trap (i.e. set n2 = Z2 = and Vb2ii, A02207 -^1221 0), which amounts to discussing two coupled 
partially-condensed systems. 

Using Eqs. (17)-(18) and (20)-(21) we can deduce the evolution of the ground state condensate mean field zo (as 
opposed to 1^0 P), which to second order reads 

dz 

— = -iujQZQ-i [Toooo - iAoooo(l + 2no)] (2no + jzop) zq + Aoooo (2no) ^0 

-2i [Toiio - 2iAoiio(l + + ni)] {ni + |zi|^) zq + 2Aoiio i^noni) zq 

-4:Xoiio {ni - no) \zif Zq (29) 

The terms in square brackets can be identified as the corresponding matrix elements of an effective many-body 
interaction introduced, strictly speaking, only over high-lying levels; this generalized effective interaction can be 
straightforwardly replaced by the usual many-body T-matrix T^jf (over all levels) , discussed in |^ , by subtracting 
from the second order expression of Eq. (29) (i.e. terms ~ ^ijji), a term corresponding to the scattering of particles 
in vacuum, so that one ends up correctly calculating the change in the effective interaction when the pair of atoms 
collides in a condensed gas (as opposed to the vacuum) |2l|j2^ , p5| , p7t . Since our analysis has been carried out in terms 
of bare particle energies LUi, such a renormalization of the terms appearing in Eq. (29) merely amounts to neglecting 
the factor of 1 in (l-|-2rto) and (l-|-no-l-rti) |3^]. We note that the same procedure should be carried out in all second 
order terms whose scattering amplitude depends on factors of the form |Ty7;;p|2:jp|2jp [{rik + l)(n-; + 1) — nf^ni] (if 
these are to be written in terms of T*^-^), and such terms appear in Eqs. (20)-(24) and (27). After 'renormalization', 
the square bracket in the first line of Eq. (29) takes the form 

f°° dr 

Toooo - jAoooo(l + 2no) = T^ooa + Tooto li™ / — e"''^ (2no) ^0^0 [Toooo] (30) 

1^0+ Jo I 

where [Tqq^q^] corresponds to the second order expression for the many-body T-matrix; a similar identification can 
be made for the square bracket of the second line of Eq. (29) (where a factor of two arises from the symmetric 
interchange of atoms and 1 corresponding to direct and exchange Hartree-Fock terms) . 

If one further ignores the last term of Eq. (29) (which corresponds diagrammatically to ignoring bubble diagrams 
with respect to the bare many-body ladder diagrams), one obtains the lowest order expression of the general equation 



= L^o^o + Tjifoo (2"o + \zo?) zo - IT^lio 



2ToYio("i + l^i|')^o-2Tc 



0110 



lim I — f 

;'?^o+ Jo I 

lim / — ( 

'?^0+ Jo « 



(„2\ rpMB^ 

y^o) 0000^0 

(2noni)roYj^zo (31) 



with the last term in each line ensuring correct scattering factors for condensate feed collisions from thermal 
atoms. This expression agrees with the equation for condensate evolution derived by one of us elsewhere Al- 
though this equation is valid when dealing with a single condensed trap level, the somewhat heuristic neglect of 
— 4Aoiio {ni — no) \zi\'^zo mentioned above suggests that it must be interpreted with some caution when the conden- 
sate spans more than one bare trap eigenstates ||3l|] . 

We now show explicitly how the above equation reduces to the coupled finite temperature nonlinear Schrodinger 
equations used for two-component condensation fl^ j. For this we must first assume that the gas is sufficiently dilute, 
so that we can further ignore in Eq. (31) the effect of the surrounding medium on binary collisions; this amounts to 
replacing Tq^^q {i = 0, 1) by the full two-body T-matrix Tq^q, while simultaneously ignoring the 'kinetic' (many-body) 
contributions corresponding to the last term of each line in Eq. (32). By further approximating T^^ by the usual 
pseudopotential UijS{r — r ) valid in 3D, where Uij is parametrized in terms of the scattering length for the collision 
of an atom in level i with an atom in level j — 0, 1) and transforming to coordinate space, we obtain 



dt 



i-^ = i?r*o + Uoo (l^oP + 2no) $0 + 2C/oi + ni) $0 (32) 
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where $i and fii denote the condensate mean field and the non-condensate density of component i, with H^^^ being 
the corresponding bare trap hamiltonian (kinetic energy plus trapping potential). By analogy, 



i-gf = i?r$i + f/ii (l^iP + 2ni) $1 + 2C/oi (l^oP + no) $i 



(33) 



Eqs. (32)- (33) correspond to the well-known finite temperatm'e Hartree-Fock equations for two-component conden- 
sation fw^ . 



2. Two-mode Condensation vs. Semi-classical Laser Theory 



In the language of quantum optics, a system in which the condensate mean field spans more than one single-particle 
eigenstate should be analogous to a photon laser in multi-mode operation. Exploring the formal connection between 
the two, with their similarities and differences, was in fact part of the motivation for this work and the main reason 
for formulating our approach in terms of bare particle eigenenergies. If such a connection were to be taken literally, 
one would expect the coupled condensate mean field equations to resemble those of the two-mode laser. Even though 
we have assumed that level 2 is not itself condensed (and may even be initially unoccupied even by thermal atoms), 
we should still also consider its presence here; this is because the existence of level 2 is an inherent property of the 
system that cannot be ignored, since it will affect the evolution of the two lowest condensed modes zq, zi via coUisional 
redistribution processes. Combining Eqs. (19)-(24) for the second order contributions to the evolution of condensate 
amplitudes in a bare single-particle basis, they are found to exhibit the general structure 



d\zo\ 
dt 



izoi^ - /3o^kor - er>oii^ii^ + e^^ikoH^^ir 



(34) 



dt 



izii^-/3-izir-0fo-izon^ii^+cfoi-onzir 



(35) 



Focusing initially only on the first three contributions of each of the above coupled equations, we note that they 
have the same form as those of the two- mode photon laser intensity equations |^^, if we identify \zi\^ with the mode 
intensity. By analogy, we thus refer to the coefficients appearing in Eqs. (34)-(35) as: the 'net gain' coefficient ai of 
each mode, the 'self-saturation' [3i of each mode and the 'cross-saturation' coefficients 9ij. In contrast to the optical 
laser where such equations arise from the polarization of the active medium, the nonlinearity in the case of the atom 
laser is intrinsic, arising from atom-atom interactions; this is absent in optical lasers, since there are no photon-photon 
interactions affecting the coherent photon field. Hence, the above coefficients a^, Pi and 9ij depend on incoherent 



populations rii and coUisional rates and ro2ii 



stand for Redistribution, 
defined by 



Self, 



The superscripts R, S, E, H used to define the above coefficients 
Exchange and Higher-level-exchange terms. The above coefficients are respectively 



RSEH 



a. 



(36) 



where 



2ro2ii [{n2 + lWi~n2{ni + lf] 
2ro2ii [non2(ni -|- 1) - (no + l)(n2 + l)7ii] 
2Aiiii [{m + - [ui + l)nf\ 
iXijji [{rii + l){nj + l)nj - {nj + l)ninj] 
■ 4:Xikki [{rii + l)(nfe + l)nfe - (rife -|- l)nink] 



(37) 



where Pq' = and 



/3f^' = (/3f)+/3f 



/3f = 2ro2ii [(no + l)(n2 -f 1) - nonj] 



/3f = Xuii [(n. + 1) 



2 - nV 



qRE _ nR I nE 



(38) 

(39) 
(40) 
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{ 




-4ro2ii [ni(n2 + 1) - (ni + 1)^2] = -6^^ 
'^^tjji {[{'"'i + + 1) - niUj] + [{rii + l)nj 



{fij + l)ni]} 



(41) 




ro2ii[(n2 + l)-«2] = -2a 



(42) 



Knowledge of their detailed form, enables us to draw important conclusions regarding the signs of the total coeffi- 
cients Q!i, (3i, 9ij and ^.y of Eqs. (34)-(35). In particular we find, just as in the photon laser, the coefficients Pi > 
always, thus giving rise to self-saturation of the mode intensity whereas coefficients ai and 9ij can be positive or 
negative (depending on the values of rii, and ro2ii). The net gain coefficient of each mode is given by Eq. (36); 
here af gives the gain coefficient due to redistribution collisions |To2iip and can be positive or negative, depending 
on the relative values of no, ni and 71.2. The remaining contributions to Eq. (36) arise from collisions of an atom 
in level i with an atom in level j or k (where i = 0, 1 and j = 0,1,2, while k ^ and leads to saturation of 

intensity growth since af , af , > always. The cross-saturation rates 9ij contain contributions from exchange and 
redistribution collisions and can be positive or negative. 

An important difference between Eqs. (34)- (35) and the corresponding ones for two-mode photon lasers is the pres- 
ence of the higher order cross-saturation terms ^ij|zop|2;i|'*, which arise solely as a result of coUisional redistribution. 
In deriving the two-mode photon laser equations, one conventionally performs a perturbative expansion of the polar- 
ization of the active medium in terms of laser intensity. For near-threshold operation and weak laser intensities, the 
third order perturbative expansion of the medium polarization is usually sufficient, and such a truncation generates 
at most terms proportional to the square of the laser intensity. Terms of third order in the laser intensity (as well as 
higher order ones) would indeed arise in two-mode photon laser theory, if one extended the perturbative treatment of 
the polarization to fifth order (or beyond). 

In the case of Bose-Einstein condensation there is no active medium to be polarized, and the nonlinearity of the 
system is due to intrinsic atom-atom interactions. In this case, the highest order jzop'^lzip" to be included in the 
two- mode equations for the 'intensities' jz^p is based on the number of single-particle operators appearing in the 
nonlinear interaction term in the hamiltonian of the system, i.e. whether one includes only two-body collisions (as 
usual, in the dilute limit na^ <^ 1), or also three-body collisions (or higher). In the usual case of two-body collisions 
as defined by Eq. (1), the resulting equations for jz^p may contain terms up to order {z*z)^; so, in general, it would 
not appear justified to ignore contributions of third order in the mode intensities. We note that the coefficients of the 
higher order cross-saturation terms appearing in Eqs. (34)-(35) have opposite signs and are given by ^10 = — 2^oi; 
where ^01 = ro2ii > 0. 

The above discussion has been given in terms of bare amplitudes |zip, and we should note that the form of these 
equations changes when the |zip refer to amplitudes of dressed modes. In particular, when the refer to amplitudes 
in modes dressed only by the condensate mean field (assuming it is physically meaningful to speak of more than a 
single condensed mode in such a dressed basis), such higher order cross-coupling contributions would not arise, thus 
yielding a direct analogy with the intensity equations for two-mode photon lasers. The justification for this is given 
in the next section, where we discuss how the second order collisional integrals become modified, upon shifting our 
single-particle eigenstates to a basis dressed by mean fields. The 'exact' analogy between two-mode condensation and 
optical lasers arising in this case might suggest that even in a basis dressed by the condensate mean field, one should, 
in principle, deal with more than one condensed modes in a trapped assembly. This analogy of a single, inherently 
multi-mode, condensate with the usual semi-classical multi-mode photon laser theory, however, appears to break down 
when shifting to a basis including higher mean field effects (i.e. those due to uncondensed atoms, anomalous averages, 
etc.) 

In this section, we have presented some similarities and differences between two-mode Bose-Einstein condensation 
and semi-classical two-mode photon laser theory. A more detailed investigation should discuss such equations in the 
presence of pumping, evaporative cooling and coherent outcoupling (in a manner analogous to single-mode atom laser 
models [^- p6| , ^ ) . The first obvious modification that would occur in this case is that the 'net gain' coefficient (s) 
of the lasing mode(s) should become positive and large, since the contribution of an irreversible evaporative cooling 
mechanism combined with the redistribution collisions should result in a large flow of particles towards such mode(s). 
At the same time, of course, the ai coefficients will acquire an additional negative contribution whose value will 
depend on the rate of outcoupling of atoms from the particular condensed level. This should enable the system 
to reach a steady state, with the coherent amplitude growth being stabilized both by the outcoupling mechanism, 
as well as by the coUisionally-induced dephasing ||37|| due to self {Xnn) and exchange {^ijji) interactions. In our 
formalism, such dephasing can be viewed as destruction of the coherent mean field amplitude |zip in favour of the 
'incoherent populations' rii (i-e. transfer of population from 'coherent' to 'incoherent'). Inherent two-body and three- 
body inelastic loss processes, which have recently been shown to be essential for reaching a steady state for an atom 
laser 133], will affect the coherent and incoherent populations of a particular mode in a different manner due to the 
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nature of our mean field decorrelation, e.g. (oJa|aJaiaiai) ~ (l^^il^ + 6nf), thus creating an essential irreversibility 
for the condensate mean field \zi\'^ to dominate over its corresponding fluctuations n^. By only outcoupling the lowest 
condensed mode of the system, one could, for example, achieve a kind of 'population inversion', in the sense that the 
largest condensate particle number accumulation might occur in some state other than the ground state of the bare 
trap. In closing this subsection, we note that the quantum nature of the multi-mode atom laser is implicit in the 
general Fokker-Planck treatment developed by Stoof [|| . 



B. Links to Other Kinetic Theories 



The expressions given earlier for the evolution of condensed and uncondensed components contain terms in the 
condensate mean field beyond order {z*z), and at first sight this appears to be in disagreement with existing kinetic 
theories P,p|,[7|,pp^ . For example, comparing our final expressions for total populations to those of Walser et al. j|] 
(in the corresponding limit of diagonal populations and no anomalous averages), we appear to obtain an excess AiV 
for the populations of levels and 2 (whereas level 1 is underestimated by twice that amount), with AA^ given by 



AiV = 2ro2ii|zi| 



A\zoY [{n2 + l)ni - 
+ koPkiP[(n2 + l 
+ |zi|4[(no + l)(n2 




(43) 



Since all terms of Eq. (43) are proportional to l^ip, it is clear that our equations will assume the usual form (e.g. as 
discussed in Walser et al. j^) in the limit of extremely weak condensation (i.e. when only the lowest trap eigenstate 
is occupied by the condensate). The natural question arising then is whether such equivalence remains beyond this 
simple limit. 

The key to understanding this apparent inconsistency is to note that our treatment has so far been given in terms of 
single-particle eigenfunctions, whose energies are bare (unshifted) ones, as if the trap were completely void. However, 
since the trap contains a large atomic medium which is condensed, the collisions will actully occur in the presence of 
the condensate, which forces the eigenenergics to vary in time due to the existing mean field potentials. To account for 
this, one can thus convert from the simple picture in terms of unshifted energies employed above, to the conventional 
one in which the effects of the mean fields are included into suitably 'renormalized' basis eigenenergies |^J^,^. In the 
next section we show explicitly that shifting our basis in a suitable manner identically reproduces the kinetic theory 
of Walser et al. j|] , thus providing an alternative derivation of the latter theory. 

To understand the relation of our theory to existing kinetic treatments, we must, in first instance, explain the 
physical origin of all terms of Eq. (43). We recall from Eq. (3) that all coUisional contributions affecting the total 

population of a particular level i arise from terms of general structure ^ 



Defining in 

the usual manner 20 2^ the quantities pn — {c\cj) and Kjk — {ckCj) and focusing initially on the first line of Eq. 
(43), we note that terms jzoPlzip arise from the elimination of p appearing in {a^a^aa) ~ [Ap{z* z) + c.c], which. 



among other contributions leads to 



dt 



Inlk [Pks{ZiZt)] 



Tklls [Pmk{ZiZt)] \ Zr + 



C.C. 



kit 



(44) 



Here the notation T^^^ stands for the corresponding two-body t-matrix multiplied by the approxi- 
mately energy-conserving integral of Eq. (10), and the above term clearly contains a contribution ^ 
{|2o2iiP^^o [("1 ~ '>^2){ziZi^)\ zi + c.c.} These latter terms would obviously not arise if the quantity p was considered 
essentially constant, as might be more appropriate, for example, for a rapidly-thermalizing 'hydrodynamic' system. 
Similarly, terms of order (i.e. second line of Eq. (43)) arise from the elimination of the pair anomalous average 
[k{z* z*) -f c.c], yielding (among other terms) 



dt 



^ ^ Tisvar (Zj Z^) I ^ ^ 



mrpq 



E 



(S) 



nilpq 



Prl 



^IrpqP" 



(ZpZq) y + C.C. 



(45) 



which leads to-- { |To2ii P(2;*zJ')(l 



no ■ 



n2){ziZi) -\- c.c.} This proves that our expressions explicitly contain anoma- 



lous condensate terms {zz) and (z*z*), in contrast to the non-condensate anomalous average k which has been 
neglected from all final expressions due to the application of the Popov approximation. 
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Finally, terms of order {z*z)^ (last contribution of Eq. (43)) can only arise from elimination of z (or z*) wherever 
it appears in the first order expressions, thus yielding 



/ J rlpq 
pql 



{Zi ZpZq 



+ z: 



T 



pql 



E 



T,, 



E 

pql 



-'-pqlsi^p^q 



pql 



C.C. 



(46) 



where the products {z* z^*'^ z) appearing immediately after each T^*'' are now the ones arising from the elimination of 
a 2; or z* from the first order expression in T{z* z* zz). This allows us to interpret the term ^ |2o2iiP|zoPl^i|* of Eq. 
(43) as arising from the terms {z*zz) generated in second order expressions, and these could be loosely interpreted 
as arising either (i) from condensate anomalous averages of the form (z^Zj), or (ii) from off-diagonal coherent terms 
(z*Zj), with the two statements being equivalent and showing that it would not be justified to include off-diagonal 
(z^Zj), while ignoring condensate anomalous averages (ziZj). 

In short, the above analysis shows clearly that all terms {z* z), (zz) and (z*z*) of second or higher order appearing 
in Eqs. (13), (17), (26) and hence (43) can only be generated by the adiabatic elimination of quanitites z, p and 
K from the corresponding first order expressions. In particular, the terms of Eq. (43) respectively arise from the 
adiabatic elimination of the quantities /0i2, Z2 and ko2- Such elimination is fully justified in a bare basis, which 
assumes that all quantities evolve in an analogous (rapid) fashion and should thus be treated consistently. However, 
tranformation to a dressed basis automatically restricts certain quantities to be slowly-evolving (which parameters are 
slowly-evolving depends explicitly on the choice of unperturbed basis dressing the single-particle eigenstates), so that 
their elimination can no longer be justified. In the next section we thus extend the above discussion to show explicitly 
that our treatment yields precisely the second order coUisional integrals of Walser et al. , upon renormalizing our 
single-particle eigenenergies to those dressed by HFB mean fields. 



1. Equivalence to Theory of Walser et al. 



To show explicitly the link of our treatment to the kinetic theory of Walser et al. Q , we now focus on the general 
expression for the normal uncondensed component in an n-level system. To establish exact analogy, in this section we 
further incorporate the uncondensed anomalous average k into our treatment, since this is explicitly present in the 
expressions of Walser et al. |^ . 

After much algebraic manipulation of the second order contributions to Eq. (26), we hence obtain the following 
result for the second order coUisional integrals arising within our bare basis 



Bare 



rp2 



Walser 



T2 



(2p. 



msPri \ ^mr^is ) 



h.c 



(47) 



The first term in the above expression has exactly the same form as the second order coUisional integrals obtained by 
Walser et al. Q (with the only difference being that in our approach the eigenenergies of Eq. (10) are now explicitly 
given by their bare values Wi), while the second term yields all the additional terms of our treatment arising when 
working within a bare-particle basis. The terms within the curly bracket appearing in the latter term denote the 
corresponding second order contributions arising from the adiabatic elimination of the quantities within the brackets, 
in the particular basis chosen. Their explicit form is obtained by the formal solution of their respective equations of 
motion given below, and such second order contributions therefore implicitly include the energy-conserving integral 
of Eq. (10). In the full hamiltonian of Eq. (1), the equations of motion for the above quantities are given by [§,0 



ifi- 



dzn 
IT 



E' — Bare < \ ^ rp r * , * 

■^nfc ^fc + ^J-nijk [Zi ZjZk -t- KjkZi 

k ij k 



2pjiZk] 



(48) 



(49) 
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(50) 



where we have ignored in the final expressions higher order correlations such as (c^cc) and (ccc) which arise in 
the above expressions (their effect will be discussed elsewhere). In the above expressions, we have further defined 
Vpq = 2 '^pkiq [Zk^i + Pik] and Ap, = Yl,ki "^piki [zkZi + Kki]- Choosing to describe our system in terms of bare trap 
eigenenergies essentially amounts to diagonalizing the bare trap hamiltonian Hq = Efj°-'^'^{alaj). In that case, the 
energies Ui appearing in the energy-conserving condition (10) correspond to the eigenstates of the harmonic trap. 

However, instead of working with a bare basis, one can choose a dressed 'unperturbed' basis {Hq + Hq) in which 
to define the renormalized single-particle eigenenergies, via 

H ={Ho+Hq) + {V -Hq) (51) 

where V corresponds to the binary collision hamiltonian defined in Eq. (1) and Hq denotes the quasiparticle hamil- 
tonian 

HQ = lY.{^P'i + ^94) + i^A'^l + ^*pqCqCp) } (52) 

pq 

with hqp = ^ {^pq^"^ + Vpq) ■ I^i this case, one is incorporating part of the coUisional evolution of Eqs. (48)-(50) into 
their 'free evolution frequencies' w^, which now become dressed to u)i. We thus find that when working in the usual 
quasiparticle basis, and upon ignoring higher order (triplet) averages, the above equations for z, p and k simplify to 
(dzn/dt) = —iujnZn, (dpji/dt) = —i {ujj — uji) pji and (dKjk/dt) = —i (Coj -\- ujk) Kjk where the uji now correspond to the 
renormalized eigenenergies in the above chosen basis. Hence (assuming as usual that initial interparticle correlations 
can be ignored), the quantities z, p and k lead to no intermediate coUisional evolution (to lowest order), which 
shows clearly that all terms of second or higher order in {z*z) or (zz) indeed vanish when dealing with single-particle 
eigenenergies dressed by the HFB mean fields; in the latter case our theory conincides with that of Walser et al (who 
actually formulated their theory in terms of unspecified number-conserving renormalized potentials as can be seen 
from Eqs. (64)-(67) of g ). 

Furthermore, this treatment brings out the implicit basis dependence of the second order coUisional integrals of 
Walser et al. In their treatment, they assume that the evolution of the system can be well parametrized by a restricted 
set of slowly-varying 'master' variables. As such, the possibility of evolution of such variables is absent from their 
theory, and hence they obtain a basis-independent formulation. Our treatment, however, indicates that their choice of 
a renormalized single-particle basis is implicit in their particular choice of the slowly-varying quantities. In a single- 
particle basis where the eigenenergies are dressed by HFB mean fields z, p and k, the coUisional integrals cannot include 
contributions due to the variations of these HFB parameters which are assumed to be static and hence our treatment 
identically coincides with that of Walser et al., defining at the same time the eigenenergies uJi as being dressed by 
HFB self-consistent potentials. This should be contrasted to the dressing caused simply by number- conserving HF 
mean fields which is implicit in Eq. (50) of their treatment. However, if one works within a bare basis, there is 
no a priori reason to assume that any quantity is slowly varying with respect to the others, and this leads to the 
additional second order contributions we have obtained in this paper (while simultaneously restricting the eigenstates 
of Eq. (10) to those of the harmonic trap). Hence, we believe our analysis further confirms the statement made in 
the last sentence of their paper, that their results should be valid when the mean-field induced energy shifts may be 
neglected during a strong coUisional event, and we thus interpret the basis correction terms of our treatment as the 
evolution of these shifts in a bare basis. Our treatment can be readily generalized to all basis-dependent shifts, by 
respectively defining other dressed basis as follows: (i) the T = Gross-Pitaevskii basis in which only the condensate 
parameters z evolve slowly and (ii) the Hartree-Fock basis in which both z and p are taken as slowly- varying, whereas 
the anomalous average k is adiabatically eliminated (even over low-lyingQ modes). Finally, in the case of an HFB 
formulation, such extra terms will only arise from triplet and higher-order averages ignored in Eqs. (48V-(50) and will 
clearly be of higher (than second) order in the interatomic potential, as also implied by Walser et al. ||]. 



^We remind the reader that the corresponding values of all quantities over high- lying modes have already been implicitly 
adiabatically eliminated, in favour of the effective two-body interaction T over low-lying modes used throughout this work. 
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V. CONCLUSIONS 



In this paper we have discussed the dynamics of the lowest-lying single-particle bare trap levels in the weak coupling 
limit, under the assumption that more than one level exhibit condensation. By comparing our equations to those 
of the two-mode photon laser, we briefly commented on the similarities and differences of such equations. We note 
that the usual discussion of two-mode condensation focuses on condensates which are either in different spin states, 
or under physical separation. Describing each of the two separate condensates by a nonlinear Schrodinger equation 
(as discussed in Sec. IV), one indeed recovers a direct analogy to the semi-classical theory of two- mode photon lasers 
(by associating the fields |<I'i(r, i)p to the laser intensity). However, since the multi-mode nature of the photon laser 
refers to different modes of the electromagnetic field, we believe that a more direct analogy can be obtained at a 
more fundamental level, namely within a single trapped condensate (see also |Q). In particular, in the case of photon 
lasers, the number of modes and the precise equations governing them depend on the atomic species, the cavity, the 
intensity and the truncation of the polarization expansion imposed. In the end, the modes into which the field will 
oscillate result from the strong coupling of the active medium to the cavity, corresponding thus, in some sense, to a 
dressed basis. In the case of a trapped Bose-Einstein condensate, the number of modes spanned by the condensate is 
determined by its size, its diluteness (i.e. whether the discussion can be essentially limited to two-body interactions at 
the particular temperatures and densities), the strength of the interactions and the trap confinement power, whereas 
the form of the equations additionally depends on the basis in which we choose to describe the system. In this case, 
the maximum factor |zip'"|zjp" which can be obtained in the equations is set by the complexity of the many-body 
interactions (i.e. two-body, three-body, etc.), whereas the factor we actually find in our final equations depends on 
the basis in which the equations are explicitly formulated. Following our renormalization discussion of Sec. V B, we 
stress that in the case of binary interactions, exact analogy with the semi-classical two-mode photon laser theory is 
obtained only within a T = Gross-Pitaevskii basis, i.e. a basis dressed only by the condensate mean field potentials 
via Hq — J2rn '^rrT'^ {o-r'^n) + (1/2) J2rsmn ^rsmnZ* ^mZn, in analogy to the abovc-mcntioned behavior of the photon 
laser. 

Changing our description slightly to explicitly include the mean-field effects on eigenenergies of our single-particle 
system, we demonstrated that our treatment is consistent with the usual description in which condensed atoms are 
described in terms of the NLSE, and the evolution of non-condensed atoms is based on the quantum Boltzmann 
equation ||^,^,^. In this case, the single-particle eigenenergies are ultimately effectively dressed by mean fields in the 
many-body T-matrix approximation, just as was found in the appr oach of Stoof |^ (with the many-body corrections 
arising from suitable inclusion of the anomalous average [p|, ^l[|25| , |2^ ). Our approach has been compared in more 
detail to the treatment of Walser et al which has been recently shown |^ to be equivalent to the Kadanoff-Baym 
Green's function formalism as applied to trapped Bose gases by Imamovic-Tomasovic and Griffin In particular, 
we have discussed how additional contributions which arise in our second order coUisional integrals upon choosing a 
simpler (than HFB) unperturbed basis, depend on the choice of this (bare or partially dressed) basis. This suggests 
an implicit assumption by Walser et al. that their eigenenergies are dressed by mean fields in the Hartree-Fock- 
Bogoliubov approximation, which is equivalent to their choice of slowly-evolving master variables. Our method thus 
yields an alternative microscopic derivation of the theory of Walser et al., based on a coupled equations of motion 
formalism. 

When applying perturbation theory only to second order in the potential, one has to justify why such a truncation is 
p hysic ally realistic. This perturbation theory is essentially a sys tem atic expansion in terms of the diluteness parameter 
Vna^. AtT — such a treatment can be justified in the limit V na^ <C 1 [Q. For a large system leading to large mean 
fields which heavily dress the single-particle eigenenergies from their bare trap values, one conventionally shifts to a 
description in terms of quasiparticles, thus employing HFB-shifted eigenener gies. In this case, the validity criterion 
of perturbation theory at finite temperature essentially becomes {kT/nUo){V na^) <^ 1 for the homogeneous system 
|^, ^j2l| , a criterion closely related to the one for the absence of critical fiuctuations occuring sufficiently close to the 
transition point p3[ |. On the other hand, the main part of this paper has been based on a treatment in terms of bare 
trap eigenenergies, since this enables a simple and direct analogy with the semi-classical equations for multi-mode 
photon laser theory. Clearly such a treatment cannot be valid for large, dense condensates, and its validity will be 
restricted to systems in which the interactions and particle numbers are so small, that the trap eigenenergies become 
only slighlty perturbed by the mean fields. A minimum (but not necessarily sufficient) criterion here is that the 
condensates are very weakly-interacting in the sense that nUo <C fiuj. By identifying the additional basis-dependent 
corrections of our treatment to the second order coUisional integrals of Walser et al, we can restate this condition as 
the requirement that the additional terms are much smaller than the first order contributions. Since the additional 
terms can be visualized as rates of change of the parameters 2, p and k, this criterion essentially reduces to the 
slow evolution of such quantities. Hence, our bare basis analysis can only be useful in the limit when all mean field 
potentials evolve very slowly, and do not heavily modify the bare trap eigenenergies. However, our explicit expressions 
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additionally include such slow evolution from a bare basis description, as opposed to existing theories in which such 
evolution is absent, applicable in the domain where the mean field energy shifts induced (on the already renormalized 
eigenenergies) during a collision can be neglected. Indeed, by re-formulating our treatment explicitly in terms of HFB 
eigenenergies, we find contributions of this type arising only in higher orders in the potential, being generated by the 
careful consideration of triplet and higher order averages, an issue which will be explicitly addressed elsewhere. 
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APPENDIX A: OFF-DIAGONAL CONTRIBUTIONS 

For completeness, we give here all off-diagonal contributions to the condensate / non-condensate populations of the 
three-level system discussed in the text. In particular, we have for the total population evolution 



(Vi - N2)pQiPio + (Nl - V2)(Coi/5io + PoiCio) 

1°^"!' \ +(iVi - No)pi2P21 + {Nl - Vo)(Cl2P21 +P12C2I 

-2ro2ii(l + 2iVi) [P02P20 + C02P20 + P02C20] 

-8ro211 {POIP12P2O + C01P12P2O + P0lCl2/320 + POIP12C20} (Al) 

/S0211 ^ _ (2) r Pw [4A'^opi2 + C12 + 2 (pioPo2 + C10P02 + P10C02)] 

°^ 1 +C10 [4A^o^Pi2 + 2Vo^Ci2 + (2pioPo2 + C10P02 + 2pioCo2)] 

_^(o) f P12 [4V2P10 + 4Vj Clo + 2 (P12P20 + C12P20 + P12C20)] 
""^^ 1 +C12 [4Vjpio + 2VjCio + (2pi2P20 + C12P20 + 2P12C20)] 

+6t?[,i^ [2VipioPi2 + 2Nl (C10P12 + P10C12) + A^i^^CioCi2] 

+e^o^ [2pioPi2 + 2 (C10P12 + P10C12) + C10C12] (A2) 

Although -RT|q2ii|2 vanishes in the limit of diagonal p, the same does not apply to Qijjl which contains a contribution 
C10C12, which will only vanish upon assuming that one of the three levels (e.g. level 2) is fully uncondensed. 
In the above expressions, in addition to the total population Ni = rii + \zi\'^ , we have defined the 'population 

terms' = i (2ni + [zip) and Nf^ = i (Sn^ -I- The population appears familiar from laser physics, as it 

contains the thermal (chaotic) contribution with a pre-factor of 2 over the corresponding condensed (ordered) term. 

Such a term will replace the full population term Ni whenever the index i is multiplied by a condensate amplitude Zi. 

By analogy, Ni — > N^'^ whenever the population term appears multiplied by two condensate mean amplitudes ZiZi 

(i.e. condensate anomalous averages of the form zz, or z*z* , but not normal averages z*z). 
In order to keep the notation general, so that it can be easily applicable to n partially condensed levels with next 

neighbour interactions, we have also defined the following rates 

"lltrj^ " ^° [4V(j=Fi)jj(j=Fi) - 2V(j=pi)(j±i)(j±i)(j^i) - V(,^i)(,zpi)(jzpi)(jzpi)] V(i=pi)(i±i)ii (A3) 

■^^li^lji = [V(jipl)ii(i;pl) + l^i(i±l)(i±l)i - Vi,iti] V(izpl)(i±l)ij (A4) 
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(i=Fl)(i±l)(i±l)(i=Fl) 



(i=Fl)(i±l)M 



(A5) 
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VkUk lim 
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For the condensed components of the two lowest states, we obtain the following off-diagonal contributions 



+4|^i|' 



,10 



^'^ooooPoiPio + (ro2ii + 2Aqiio) P12P21 + 2A0QOQP02P20 
+2(27o%^-7o%^)PioPi2 

[AoiloPOlPlO + ro21lPl2P21 + 270110P10P12] 
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,,0211 



7o"o'o^o')pi2] 
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-4/301 [ro2ii('^i - «2) - 2Aoiioni] 

-4pi2 [-7o°o'o'o'no + (7o°?!o' + H'Al - H¥2l) ni - 

-4 [(A0211 + 2Ag?fg) P02P21 + 70000V10/902] 

" 270000 C10P02 



. (470211 



'1221 



-72°l2^2^)n2] 



(A9) 



Xii(y2) = -4|0i|2| 



(ro2ii + 2A?}i;) poiPio + (r, 



^ ,^ J y iiii, . V,., , 0211 + 2A1111) P12P21 + (-ro211 + 2A0110) P02P2O 

1 1 -1-9 ('9-^,0211 I 9„,0211 „,0211^ „ ^ 

L +2 (,270110 + 271221 - 7iiii j P10P12 

+4|2;o|^ [AoilO/OOlPlO + ro211/9l2P21 + 27oiJoPlOPl2 
f nl.. |2 r/o-n , ,<\01in\ „ , ;'r-0211 , 



} 



+C10 < 
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+4/901 

-4pi2 
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Although it is cxtcrcmcly hard to discuss the physical implication of each separate contribution appearing in the 
above equations, at this point we would like to comment briefly on a carefully selected subset of the above equations, 
namely 



71 1 2 

= 8\zi\'^ {AoiioPoiPlO + ro21lPl2P21 + 27^iJoPlOPl2] 



(All) 



d\zi\ 
dt 



■■ 8|2:o|^ {AoiloPOlPlO + ro211/9l2P21 + 270110PIOP12] 



(A12) 



Such terms indicate clearly a mechanism of 'coherent population transfer', that is growth of condensation in one (bare) 
level due to the existence of condensation in another (bare) level, which occurs via off-diagonal incoherent couplings 
Pij. By further ignoring, for simplicity, the redistributional processes (V0211 — * 0) we see a sort of Rabi-like oscillation 
between the two condensates, via the process 



d\i 



dt 



[PijPji] 



(A13) 
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Such processes do not lead to changes in the total trap level populations. On the contrary, due to the co-existence 
of coherent and incoherent atoms within the same trap level, this term can be interpreted as follows: If one level 
initially contains a non-vanishing coherent amplitude, this will be gradually transferred to the other level (even if 
the other level is initially fully incoherent) without an associated change in the total population of each level (since 
the incoherent population of each level adjusts accordingly to keep Ni fixed). Hence, even if we assume that the 
condensate initially resides only on the bottom trap level, the above contribution generated by nonlinear coupling 
interactions will tend to give rise to a coherent mean field amplitude in level 1 (and vice versa). Whether this term 
actually becomes important at all (and in what limits this may be so) depends on how heavily this contribution 
is overshadowed by all other existing terms in the equations of motion for coherent evolution. Although we would 
not expect this term to play a significant role in the full population dynamics, this is something which should be 
confirmed by direct numerical simulation of the full equations given in this paper. We conclude the discussion by 
noting that such processes may become important in the case of coupled 2-species condensation or in the 

creation of non-ground state trapped condensates p5| , whereby the population between such levels is controlled by 
the application of external fields. 

The off-diagonal contributions to the uncondcnsed dynamics of the two lowest levels are given by 
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and 
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APPENDIX B: DERIVATION OF CONDENSATE EVOLUTION 



The second order contribution to |z,;|2 of Eq. (17) can be found by adiabatically eliminating the quantity z*(aja„iar 



via 
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Hence, the quantity F(^i^\i^i^^ corresponding to F of Eq. (||), now becomes 
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APPENDIX C: SECOND ORDER COLLISIONAL INTEGRALS 

For completeness, we give here the second order coUisional integrals of Walser et al. [|| as generated by our approach. 



Wals 



As explained in the text 
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Here Alu = (lui + ut 



X^^j — I ill Q-'i-('^i+'^t-'^p-i^q){t-t ) _ 



= tt6{Auj) - iP 



(C2) 



and the uji denotes the eigenenergy of level i in the particular basis chosen for the 



analysis of the system. In our original formulation, these correspond to bare trap eigenenergies, whereas in order to 
establish exact analogy with the coUisional integrals of Walser et al. , these should be replaced by their corresponding 
values dressed by normal and anomalous HFB mean fields. 

In the original formulation of Walser et al., their basis is left unspecified, and their resulting expressions are assumed 
to be valid for any single-particle basis dressed by, at most, number- conserving mean fields (i.e. no anomalous 
averages), a conclusion not supported by our analysis. 
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